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ABSTRACT 


Velocity  and  temperature  fluctuations  from  multi-level 
measurements  in  the  marine  boundary  layer  were  examined  to 
determine  the  validity  of  present  formulations  and  prediction 
techniques,  Spectral  distributions  were  analyzed  using  an 
analog  analyzer  as  opposed  to  the  digital  computer  methods 
previously  utilized  at  the  Naval  Postgraduate  School. 

In  general,  the  height  variations  of  the  dissipation 
rate  of  turbulent  kinetic  energy,  momentum  flux,  and  the 
temperature  structure  parameter  observed  in  this  study  agree 
with  the  overland  expressions  developed  by  Wyngaard,   Most 
deviations  from  expected  results  were  attributed  to  apparent 
anomalies  in  spectral  results. 

It  has  been  shown  that  the  analog  spectral  analysis 
performed  in  this  study  produced  information  which  is  as 
reliable  as  that  from  digital  methods, 
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I,   INTRODUCTION 

This  is  a  study  of  velocity  and  temperature  fluctuations 
in  the  marine  boundary  layer.   The  data  examined  were  obtained 
from  multi-level  measurements  with  tungsten  wire  velocity 
probes  and  platinum  wire  temperature  probes.   The  measure- 
ments were  made  from  the  Naval  Postgraduate  School  research 
vessel,  ACANIA,   Mean  wind  and  temperature  measurements  are 
also  utilized  to  complement  and  analyze  the  fluctuation  data. 
A  separate  investigation  on  the  mean  data  has  been  completed 
by  Smedley  (1975). 

The  purpose  of  this  study  is  to  examine  properties  of 
small  scale  temperature  and  velocity  fluctuations  in  the 
marine  environment,  on  the  basis  of  spectral  distributions. 
A  significant  aspect  of  the  study  is  that  spectral  distribu- 
tions were  obtained  utilizing  an  analog  spectrum  analyzer, 
rather  than  digital  computer  methods.   The  latter  has  been 
the  method  used  previously  in  atmospheric  turbulence  experi- 
ments at  the  Naval  Postgraduate  School.   The  study  includes 
descriptions  of  the  variation  with  height  of  (1)  z    (the  viscous 
dissipation  of  turbulent  kinetic  energy)  and  (2)  ci  (the 
temperature  structure  parameter). 
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II.   THEORETICAL  CONSIDERATIONS 

A .  BACKGROUND 

This  study  is  one  of  a  continuing  investigation  of 
multiple-level  overwater  observations  from  a  ship.   The  ship 
is  the  R/V  ACANIA,  shown  in  Figure  1,  operated  by  the  Depart- 
ment of  Oceanography,  Naval  Postgraduate  School.   Because 
these  were  shipboard  measurements,  only  high  frequency  small 
scale  fluctuations  could  be  measured,   The  following  paragraphs 
will  review  the  theoretical  considerations  used  to  interpret 
the  small  fluctuation  data,   Parallel  reviews  of  this  theory 
have  been  provided  by  Johnston  (1974)  and  Bone  (1974).   The 
latter  involved  initial  examinations  of  data  obtained  with 
the  same  instrumentation  used  in  this  study, 

B.  TURBULENT  KINETIC  ENERGY  BALANCE  EXPRESSION 

In  a  simplified  form,  the  turbulent  kinetic  energy 
balance  expression  for  neutral  conditions  can  be  written 


-u'w'   9u  z- x 

3i  =  E  (1> 

where  u  represents  the  mean  wind  at  height  z  and  the  term 

y-  is  the  wind  shear,   Equation  (1)  states  that  mechanical 

production  of  turbulent  kinetic  energy,  -u'w'  y- ,  equals  the 
molecular  dissipation  of  turbulent  kenetic  energy,  ?• 
Equation  (1)  is  often  written  as 

u2  A"  =  £  (2) 
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Figure  1.   The  R/V  ACANIA 
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where  u£  =  -u'w* ,  and  -u'w'  represents  the  vertical  transfer 
of  horizontal  momentum  by  turbulent  fluctuations. 

For  neutral  conditions  in  a  constant  stress  boundary  layer, 

the  following  expression  is  valid  for  the  wind  gradient 

~-   u 

—  =  —  (3) 

8z    Kz  (6) 

where  K  is  Von  Karman ' s  constant.   Combining  Equations  (2) 

and  (3)  yields  an  expression  to  estimate  u   on  the  basis  of 

measurements  of  e  at  a  known  height  z. 

u£  =  eKz  (4) 

Equation  (4)  can  be  rewritten 

Ln  £  =  -Ln  z  +  Ln  =-  (5) 

which  is  a  form  used  in  later  analysis  of  the  variation  of 
e  with  height. 

For  non-neutral  conditions  the  turbulent  kinetic  energy 
balance  expression  (Garratt  1972)  is 


fon  - 


For  an  unstable  case  (R.  <  0) , 

*  (€)  -  (1  -  16  £)"*  (7) 

m 

and  £  =  R.  . 

For   the   stable   case    (R.    >   0), 

l    .  ' 

<frmU)   =    1   +   5    5  (8) 

and 

R. 

5  « 


1    -    5   R. 

l 
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where  R.  is  the  Richardson  number. 
1  ' 

R,  =  g/T  (3T/3z)  /  (au/iB)« 
Equation  (5)  can  be  rewritten  as 

u3 
Ln  c  =  -Ln  (z/(vo  _5))+Ln  _f  <9> 

A  plot  of  Ln  e  versus  the  first  term  on  the  right  hand 
side  of  Equations  (5)  or  (9)  would  yield  a  straight  line  with 
-1  slope.   Figure  2  illustrates  the  possible  effects  of 
instability  on  the  predicted  plot  of  Ln  z  versus  Ln  e,  where 
Ln  <j>  (R.)  is  an  empirical  function  of  the  Richardson  number, 
Deviations  from  this  relationship  could  be  interpreted  as 
either  due  to  deviations  from  the  constant  flux  assumption 
arising  from  the  influence  of  wind-wave  coupling  or  non- 
neutral  stratifications. 

From  the  above  discussion,  e  can  be  used  to  estimate  u^ 
values.,  u^  is  an  important  parameter  in  boundary  layer  studies 
since  it  defines  the  intensity  of  turbulent  transfer  of 
momentum  toward  the  surface.   It  is  assumed  to  be  constant 
with  height  in  the  first  30  meters. 

C.   SPECTRAL  SIMILARITY  THEORY 

Kolmogorov's  second  hypothesis  is  the  basis  for  estimating 
£  and  C*  (the  temperature  structure  function  parameter) 
from  spectra  of  velocity  and  temperature  fluctuations.   It 
is  based  on  similarity  expressions  for  the  inertial  subrange, 
which  is  the  wave  number  region  where  energy  is  transferred 
without  dissipation  from  lower  to  higher  wave  numbers. 
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Figure  2,   An  Illustration  of  the  Effect  of  Stability  on 
Variation  of  e  with  Height, 
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The  following  relationships  were  proposed  by  Kolmogorov 
for  spectral  densities  in  the  inertial  subrange, 

4>(k)  =  Cj  e  2/3  k"5/3  (10) 

4>T(k)  =  c2  C*  k"5/3  (11) 

where  $(k)  and  <f>T(k)  are  the  one  dimenensional  spectra  for 
velocity  and  temperature  respectively  and  k  is  the  streamwise 
component  of  wave  number  as  defined  below. 

e  has  already  been  discussed  in  it's  relationship  to  u+. 
Another  turbulence  parameter  of  interest  in  this  study  is  the 
temperature-structure  parameter,  C*   which  appears  in 
equation  (11).   It  is  important  because  it  is  related  to  the 
refractive-index-structure  parameter,  C*   by  the  expression 

ca  .  /79  x  10-'P y    ca  (12) 


,2    /79 


.2 


T 


C^  is  a  primary  parameter  in  the  characterization  of  optical 
propagation  in  the  environment. 

Height  variations  of  C£  are  dependent  on  stability,  as 
was  the  case  for  e.   Predicted  height  variations  for  various 
stability  conditions  are  shown  in  Figure  3.   Wyngaard,  et .  al 

—  2/3 

(1971),  noted  that  C3,  decreased  with  height  as  Z*"  '"  for 
stable  conditions  and  as  Z~h'       for  unstable  conditions. 

Since  velocity  and  temperature  fluctuations  are  measured 
at  a  fixed  point  in  the  flow,  the  resultant  spectra  are 
realized  at  a  "temporal"  frequency,  f,   To  obtain  an  e  and 
C™,  the  time  (f)  and. space  (k)  scales  must  be  related  by 
Taylor's  (1938)  "frozen  turbulence"  hypothesis,  i.e., 
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Figure  3.   An  illustration  of  the  effect  of  stability  on 
the  variation  of  Ci,  with  height. 
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k  =  27Tf 
u 

where  u  is  the  mean  wind  speed  at  the  measurement  level.   The 

term  "frozen"  implies  the  turbulence  pattern  remains  unchanged 

as  it  sweeps  past  the  probe. 

Equations  (11)  and  (12)  can  also  be  written  as 

f(J>(f)  =  k<J>(k)  =  c  e2/3k~2/3  (13) 

f<f>T(f)  =  k<f>T(k)  =  c2c^k-2/3  (14) 

where  <t>(f)    and  <{>T(f)  are  spectral  density  values  with  units 
of  (m/sec)2/Hz  and  C2/Hz  respectively, 

Final  forms  of  Equations  (13)  and  (14),  in  terms  of  4>(f) 
and  <f>  (f),  which  are  measured,  are 

/  1   2  /  3  \  3  /  2 

E  -  {  ^ f<Kf)  (15) 


^ 


C2  =     k^/3   fA   (f )  (16) 

c 

2 

where,  empirically,  c l~  0.5  and  c2  =  0.25. 

From  Equations  (15)  and  (16),  with  measured  f,  <j>(f),  ^^(f), 
and  u,  E and  C2  can  be  determined  at  each  level  of  interest. 
These  were  the  final  forms  used  to  determine  the  variation 
with  height  of  e  and  C2  values. 
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III.   DATA  COLLECTION  AND  SENSOR  CALIBRATION 

A,  THE  PLATFORM 

Observations  were  obtained  aboard  the  R/V  ACANIA  anchored 
off  Monterey,  California  in  Monterey  Bay,   Measurements  were 
made  at  four  levels  on  two  masts  spatially  mounted  on  the 
forward  deck  of  the  ship.   Figure  4  illustrates  the  sensor 
locations  on  these  masts.   Figure  5  shows  the  vane  and  probe 
arrangement,   The  vane  maintained  a  one  dimensional  profile 
of  the  wind  during  measurements, 

B.  INSTRUMENTATION 

The  velocity  fluctuation  sensors  were  TSI  model  1210 
probes  with  sliding  support  shields  and  tungsten  wires.   The 
shields  permitted  isolation  of  the  sensing  region  for  deter- 
mination of  the  undisturbed  velocity,  V  ,  before  and  after 
each  experiment.   The  platinum  coated  tungsten  wire  was  small 
enough  to  resolve  the  viscous  dissipation  scale  without  any 
need  to  apply  wire  length  corrections.   The  wire  was  fitted 
with  plated  ends  for  isolating  the  sensing  area  and  thereby 
minimizing  flow  disturbance. 

Electronics  associated  with  each  probe  were  a  TSI  Model 
1054B  Linearized  Anemometer  ,   a  TSI  Model  1056  Variable  Decade 
Module  and  a  Model  1057  Signal  Conditioner.   The  signal 
conditioner  was  only  used  during  the  March  1974  observations. 
The  anemometer  had  a  linear  frequency  response  from  DC  to  10 
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KHz  and  the  variable  decade  module  operated  with  a  0-60  ohm 
range.   The  signal  conditioner  had  a  linear  frequency  response 
from  DC  to  400  KHz  and  permitted  voltage  suppression  in  one 
volt  steps  from  0-29  volts. 

Sensor  placement  required  exceptionally  long  cables,  but 
any  decrease  in  system  frequency  response  had  little  effect 
in  the  frequency  band  of  interest, 

In  general  the  performance  of  the  hot-wire  system  was 
excellent  throughout  the  experiments. 

Temperature  fluctuations  were  also  measured  with  a  TSI 
Model  1210  configuration  using  platinum  wire  sensors. 
Electronics  used  were  the  GTE  Sylvania  Model  140  thermosonde 
temperature  modules,  which  can  be  used  to  make  both  single 
point  temperature  fluctuation  and  AT  measurements.   The 
temperature  module  is  designed  for  measurements  of  temperature 
variations  down  to  0,001°C  at  frequencies  up  to  1  KHz. 

Mean  measurement  systems  are  described  in  Smedley's 
(1975)  study. 

C.   SENSOR  CALIBRATION 

1.   Hot  Wires  Calibration 

In-situ  calibrations  of  the  velocity  sensors  were 
accomplished  several  times  during  each  experiment,   In  this 
procedure,  recordings  v/ere  made  of  both  the  cup  anemometer 
RPM  count  and  the  corresponding  hot  wire  voltage  output.   The 
sensor  wind  speed  and  voltage  relationship  is  given  by 

v2  =  a(u)*  +  b  (17) 
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where  v  is  the  voltage  and  u  is  mean  wind  speed  for  a  partic- 
ular level.   The  constants  a  and  b  are  calibration  curve  slope 
and  intercept  respectively  obtained  from  the  in-situ  calibra- 
tion.  Figure  6  is  a  sample  calibration  curve  constructed 
from  in-situ  values  in  one  period. 

To  convert  the  power  spectral  density  levels  to 
velocity  units  required  a  calibration  factor  such  that 

c  v»  =  uf  (18) 

where  c  is  the  calibration  factor  in  — '  , .  ,  v'  is  the 

volt  ' 

voltage  flucuation  and  u'  is  the  velocity  fluctuation, 
Differentiating  Equation  (17),  we  get 

(ivj£)i)  V  -„•  (19) 

or 

c  =  iz£l*  (20) 

a 

To  obtain  the  calibration  factor,  c,  Equation  (17) 
was  initially  solved  for  v  using  the  appropriate  sensor 
level's  mean  wind  speed  u  and  also  the  particular  sensor's 
calibration  curve  slope  and  intercept. 
2.   Platinum  Wires  Calibration 

Wire  resistance  is  related  to  temperature  by  the 
equation 

R  =  Ro(l  +  ci(T  -  To))  (21) 

where  R  =  resistance,  T  =  calibration  temperature,  To  =  0°C, 
Ro  =  resistance  at  To,  and  a  =  ,0039/°C  (temperature 
coefficient  of  resistance  for  platinum).   Each  wire's 
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Figure   6.      Hot   wire   calibration   plot. 
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resistance  was  measured  in  the  laboratory  at  a  temperature  T, 
to  determine  Ro, 

We  can  relate  temperature  fluctuations  to  wire 
resistance  fluctuations  by  differentiating  Equation  (21), 

dR  =  a  Ro  dT 
or 

R'  =  a  Ro  T'  (22) 

Recorded  voltage  fluctuations  can  then  be  related  to 
resistance  changes  by  the  transfer  function,  H,  in  the 
expression 

v»  =  H  R'  (23) 

where  H  was  determined  for  the  four  bridges  by  varying  the 
resistance  of  a  TSI  1056  Variable  Decade  Probe  Resistance 
and  recording  the  corresponding  voltage  changes,   Figure  7 
is  a  calibration  curve  of  voltage  versus  resistance,  where 
the  slope  determines  the  transfer  function,  H. 

Finally,  to  relate  recorded  voltage  fluctuations  to 
temperature  fluctuations  substitute  Equation  (23)  into 
Equation  (22),  giving 

T»  =  I 1 ]  v  (24) 

\  H  a  Ro  j  v   ' 

or 

T'  =  Cu   v»  (25) 

where  C„  is  the  temperature  calibration  factor  in   C/volt. 
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Figure  7.   Platinum  wire  transfer  function  plot. 
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D.   SIGNAL  PROCESSING 

The  data  from  the  four  levels  was  recorded  simultaneously, 
in  analog  form,  on  a  14  channel  tape  recorder. 

The  signals  were  continuously  monitored  with  an 
oscilloscope  and  an  8  channel  strip  chart  recorder  to  insure 
the  quality  of  the  recorded  signal  and  to  detect  equipment 
failures. 
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IV,   DATA  PROCESSING  AND  ANALYSIS 

A,  PRELIMINARY  ANALYSIS 

Selection  of  data  periods  was  based  on  signal  evaluations 
of  strip  chart  records,  mentioned  in  section  III.D,  and  on 
evaluations  of  experimental  logs  of  equipment  status  and 
ship  positioning.   Primary  consideration  was  given  initially 
to  obtaining  18  minute  records  with  temperature  and  velocity 
signals  at  all  levels,   This  was  later  changed  to  21  minute 
records  due  to  characteristics  of  the  anolog  spectrum  analyzer, 

B .  PROCEDURES 

Several  procedures  were  adopted  to  best  handle  the 
turbulence  data.   These  are  outlined  in  Appendix  A  in  detail. 
All  spectral  analysis  was  performed  by  processing  the  recorded 
analog  signal  on  a  Federal  Scientific  Model  UA-500-1  Spectrum 
Analyzer-Averager  shown  in  Figure  8. 

The  averager  output  was  a  logarithmic  plot  of  the  power 
spectral  density  in  Volts2/Hz,  versus  frequency  f  in  Hz.   The 
plotter  was  calibrated  so  that  all  plots  would  concide  with 
Bone's  (1974)  digital  plots.   Thus,  a  preliminary  visual 
comparison  was  made  to  verify  the  spectra  prior  to  computations. 
Figure  9  illustrates  a  sample  of  the  plot  comparison  for  a 
data  period  common  to  Bone's  (1974)  study  and  this  study. 
Figure  10  shows  a  typically  good  graphical  plot  of  spectrum 
versus  frequency  used  in  this  study.   Figure  11  is  an  example 
of  a  poor  spectral  plot. 
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Figure  9,   Comparison  of  digital  and  spectrum  analyzer 
velocity  spectra  for  1855-1913L,  18  Jan.  74 
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Figure  10,   Typical  good  analyzer  spectral  plot, 
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Figure  11,   Typical  poor  analyzer  spectral  plot, 
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C,   SPECTRAL  INTERPRETATION 

Scaling  of  the  spectral  plot  to  relate  an  RMS  input 
voltage  to  a  power  per  unit  of  frequency  is  discussed  in 
Appendix  B. 

The  single  point  method  as  previously  discussed  by 
Johnston  (1974)  and  Bone  (1974),  was  used  to  obtain  e  and 
Ci,  values  from  the  velocity  and  temperature  spectra.   Using 
many  spectral  points  in  the  inertial  subrange  (-5/3  slope 
region),  an  intercept  at  an  arbitrary  frequency  f  can  be 
obtained,   This  intercept  value  yields  the  power  spectral 
density  (PSD),  after  correction. 

1.   Determination  of  c    and  u 


The  measured  PSD  value  was  converted  to  a  spectral 
density  in  engineering  units  by  the  equation 

(J)(f)  =  c2  •  PSD 

=  (cm/sec/volt)2  •  volt2/Hz  -  (cm/sec)2/Hz  (26) 
where  c  is  the  hot-wire  calibration  factor.   The  solution  to 
Equation  (15)  yielded  a  value  for  z   at  each  level.   Then, 
Equation  (4)  and  (6)  were  solved  for  a  u+  value  and  a 
corrected  u^  value  respectively.   Corrected  u^  values  were 
not  valid  for  Richardson  numbers  >  0.2. 
2.   Determination  of  C2 

The  following  equation  was  used  to  convert  PSD 
values  from  measured  spectra  to  dimensional  values 

4>T(f)  =  C£  •  PSD 

"  <°C/volt>2'  Volt2/Hz  =  °C2/Hz  (27) 
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where  C„  is  the  calibration  factor  for  the  temperature 
n 

system.   Solving  Equation  (16)  for  C*  yielded  a  value  for 
each  level  Z,   These  values  were  then  applied  to  Equation 
(12)  to  obtain  C^  for  possible  comparison  with  optically 
derived  Cil  values, 
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V,       RESULTS 

A,  DATA  PERIODS  EXAMINED 

Fluctuating  velocity,  u',  and  temperature,  T',  data  were 
studied  for  the  nine  periods  shown  in  Table  I ,   Four  levels 
of  both  velocity  and  temperature  were  desired,  but  temperature 
signals  were  not  always  available  for  all  levels  due  to 
various  component  failures,-   Table  I  summarizes  the  periods 
examined  and  what  data  was  available  for  each  period, 

B.  VELOCITY  SPECTRAL  RESULTS 

For  seven  of  nine  periods  studied,  four  levels  of 
velocity,  u' ,  were  available.   Equation  (15)  was  used  to 
calculate  z    from  values  at  each  level  from  selected  frequencies, 
f,  spectral  densities,  $(f),    and  wind  speed,  u,   Each  computed 
value  e  was  then  used  in  Equations  (4)  and  (6)  to  determine 
a  value  for  the  friction  velocity,  u^, 

Table  II  summarizes  the  £  and  u^  results  along  with 
stability  corrected  values  of  u+,  profile  u+  estimates 
obtained  by  Smedley  (1975),  and  Richardson  numbers,  R . .   Most 
spectra  exhibited  the  -5/3  slopes  for  the  inertial  subrange 
region  which  validates  the  method  used  to  determine  £  from 
the  spectral  densities.   Except  for  one  value  on  March  27, 
u^  values  are  all  less  than  0.5  m/sec-   and  most  are  in  the 
expected  range  of  0.1  -  0.5.   Most  periods  of  u^  seem  to 
reflect  the  constant  flux  theory  accurately,  especially  if 
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TABLE  I 


DATA  PERIODS  IN  THIS  STUDY 


PERIOD 


DATE 


TIME 


LEVEL 


PARAMETER 


8 


27  Mar, 

74 

1841-1902L 

1 

T1 

2 

u' 

,  T' 

3 

u  ■ 

,  T» 

• 

4 

u' 

27  Mar, 

74 

1902-1923L 

1 

u' 

,  T' 

2 

u' 

,  T' 

3 

u' 

,  T' 

4 

u ' 

27  Mar, 

74 

1928-1945L 

1 

u ' 

,  T' 

2 

u' 

,  T' 

3 

u ' 

,  T* 

4 

u' 

,  T' 

27  Mar, 

74 

2117-2138L 

1 

u' 

,  T' 

2 

u' 

,  T* 

3 

- 

4 

u' 

,  T' 

28  Mar. 

74 

1902^1923 

1 
2 

u ' 
u' 

,  T' 

3 

u' 

,  T' 

4 

u1 

,  T' 

28  Mar. 

74 

2059-2117L 

1 

u' 

,  T' 

2 

u ' 

,  T' 

3 

u' 

4 

u ' 

,  T' 

18  Sept 

.  74 

1643-1701L 

1 

u' 

,  T' 

2 

u' 

,  T' 

3 

u' 

4 

u' 

,  T' 

19  Sept 

.  74 

1735-1753L 

1 

u* 

,  T' 

2 

u' 

,  T' 

3 

u' 

4 

u? 

,  T' 

19  Sept 

.  74 

1900-1918L 

1 

u' 

,    T' 

2 

u' 

,  T' 

3 

u' 

4 

u' 

,  T' 
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TABLE  II 


e  and  u^  Results 


HEIGHT 
(m) 

£ 

(m2-secT3  ) 

1           2 

u*     _     u* 
(m-sec.1)    (corrected) 

u* 

( profile) 

R. 

l 

Period 

1 

27 

March 

74 

1841-1902L 

6,6 

7,6 

13.8 

,0061 
.0069 
.0069 

,25 
.27 
.33 

* 

* 
* 

.15 

.53 
.62 
.43 

Period 

2 

27 

.March 

74 

1902-1923L 

4.1 

6.6 

7.6 

13.8 

.021 
.012 
.019 
.013 

.32 

.31 
.38 
.41 

.29 
.19 
.28 
* 

.25 

,06 
.14 
.16 
.38 

Period 

3 

27 

March 

74 

1928-1945L 

4.1 

6.6 

7.6 

13.8 

.240 
.039 
,018 
.021 

.72 

.46 
.37 
.48 

.65 
.34 
,25 

* 

.26 

,06 
,13 
.15 
.35 

Period 

4 

27 

March 

74 

2117-2138L 

4.1 

6.6 

13.8 

.044 
,025 
.016 

.41 
,40 
.44 

.34 
.36 
,25 

.39 

,09 
.06 
.17 

Period 

5 

28 

March 

74 

1902-1923L 

4.1 

6,6 

7.6 

13.8 

.039 
.0030 
.0030 
.0013 

.40 
.20 
.20 
.19 

.38 
.18 
.19 
.14 

,30 

.02 
.05 
.05 
.13 

NOTE 

1 

computed  from  equation  15 


computed  from  equation  6 
*  no  data 
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TABLE    II    (continued) 


HEIGHT 

E 

u* 

u* 

u* 

R. 

l 

(corr . ) 

(profile) 

Period 

6 

28 

March  74 

2059-2117L 

4.1 

,020 

.32 

.27 

.17 

.09 

6.6 

.013 

.32 

,14 

.19 

7.6 

.010 

,31 

* 

.22 

13.8 

.011 

,38 

* 

.52 

Period 

7 

18 

September 

74 

1643-1701L 

4.1 

.0089 

.24 

.19 

.10 

.12 

6.6 

.0015 

.16 

* 

.22 

7.6 

.00049 

.11 

* 

.26 

13.8 

.019 

,46 

* 

.62 

Period 

8 

19 

September 

74 

1735-1753L 

4.1 

.0099 

,25 

.26 

.16 

-.11 

6.6 

.0063 

.25 

.26 

-.23 

7.6 

.0017 

.17 

.17 

-.27 

13.8 

.0039 

,27 

.26 

-.64 

Period 

9 

19 

September 

74 

1900-1918L 

4.1 

.0077 

,23 

.24 

.25 

-.04 

6.6 

.0076 

.27 

,28 

-.08 

7.6 

.0013 

.15 

.16 

-.09 

13.8 

.045 

,62 

.64 

-.22 
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obviously  bad  levels  are  discarded  in  the  periods,   Those 
results  which  prevent  nearly  the  same  u+  flux  for  all 
levels  in  a  period,  were  readily  attributed  to  a  departure 
from  the  -5/3  spectral  slope,   A  comparison  of  this  study's 
u+  and  stability  corrected  u^  values  with  Smedley's  regression 
u^  values,  provides  close  agreement  for  at  least  two  levels 
for  each  period,   The  comparatively  high  u+  values  for  level 
4  in  the  stable  cases  has  not  been  explained,  but  it  also 
occurred  in  one  of  the  two  unstable  cases. 

Figure  12  is  a  typical  plot  of  height  versus  e  for  a 
stable  case,   It  shows  close  agreement  with  the  -1  slope 
expected  for  the  near  neutral  case,  according  to  Equation  (5), 
Figure  13  illustrates  the  unstable  case.   Here  again  we  have 
agreement  with  the  neutral  case  -1  slope,  except  for  level 
3.   Poor  results  for  level  3,  in  this  case,  can  be  related 
to  a  poor  spectral  plot  with  insufficient  amplitude  and  a 
slight  departure  from  a  -5/3  slope.   This  deviation  is  not 
expected  since  level  2  and  3  are  within  one  meter  in  height. 

Composites  of  height  versus  e  for  all  periods  are  shown 
in  Figure  14  for  stable  and  unstable  cases.  Obviously  bad 
data  was  not  included  in  these  composites.  Both  the  stable 
and  unstable  cases  show  close  agreement  to  the  near  neutral 
prediction.  This  differs  from  the  previous  results  of 
Bone  (1974),  in  that  his  slopes  deviated  from  the  -1  value, 
as  expected  for  non-neutral  conditions. 
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C,   TEMPERATURE  SPECTRAL  RESULTS 

One  period  had  four  levels  of  temperature,  T' ,  data, 
while  eight  periods  only  had  three  levels  of  T'  data.   Ci, 
values  were  computed  by  applying  Equation  (16)  to  values  at 
each  level  for  selected  frequencies,  f,  spectral  densities, 
<}>T(f),  and  wind  speeds,  u,   Equation  (12)  was  then  applied 
to  the  resultant  C,i  to  get  CA,   The  C*  values  will  be  used 
by  investigations  into  coincident  optical  propogation 
experiments , 

Table  III  is  a  summary  of  the  Ci,  and  C^  results  for  this 
study.   It  also  includes  Smedley's  (1975)  profile  values  of 
Ci,  for  a  comparison.   The  height  variations  of  Ci,  in  the 
stable  periods  were  generally  consistent  with  predictions. 
During  unstable  pexiodt>  there  appealed  to  be  some  lack  of 
consistency  between  levels.   Here  again,  most  inconsistencies 
were  traced  back  to  a  poor  spectra,   Ci,  values  which  are 
within  an  order  of  magnitude  of  Smedley's  results  are  consid- 
ered a  fair  correlation.   Such  correlation  does  not  occur  for 
any  level  in  periods  1,  4,  and  7, 

A  typical  plot  of  height  versus  Ci,  for  the  stable  case 
is  shown  in  Figure  15.   It  agrees  very  well  with  the  expected 
-3/2  slope,   The  unstable  case  is  shown  in  Figure  16,   The 
expected  slope  of  -3/4  for  height  versus  Ci,  agrees  reasonably 
well  with  two  of  the  three  data  points.   The  4  meter  level 
point  can  be  attributed  to  a  spectra  with  a  slope  greater 
than  -5/3  along  with  comparatively  small  amplitudes. 
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TABLE  III 


C2  and  C2T  Results 


'N 


HEIGHT 
(m) 

c21 

UT 
(°C2  m' 

■«/•) 

c22 

CT 
(profil 

e) 

C2 
N 
/  -2/3  -,  rt"1  1  x 

(m   '  xlO    ) 

R. 

l 

Period 

1 

27 

March 

74 

1841. 

-1902L 

4,1 
6.6 
7.6 

.0017 
.0034 
,0034 

* 
* 
* 

.16 
.32 
.33 

.25 
.53 
.62 

Period 

2 

27 

March 

74 

1902- 

-1923L 

4,1 
6.6 
7,6 

.0021 
.0022 
.0043 

,043 

,0061 

.0032 

,20 
.21 
.42 

.06 
.14 
.16 

Period 

3 

27 

March 

74 

1928- 

-1945L 

4.1 

6.6 

7.6 

13.8 

.0057 
.0050 
.0062 
.0015 

.042 
.0065 
,0036 
* 

.55 
.48 
.60 
,14 

.06 
.13 
.15 
,35 

Period 

4 

27 

March 

74 

2117- 

-2138L 

4,1 

6.6 

13.8 

.0036 
.0032 
.0018 

,072 
,020 
.00089 

.35 
.31 
.17 

.09 
.06 
.17 

Period 

5 

28 

March 

74 

1902- 

-1923L 

4.1 

7,6 

13.8 

.017 
.026 
.0018 

,016 

.0040 

,00055 

1.7 
2,6 

.18 

.02 
.05 
.13 

value  computed  from  Equation  (16) 
values  from  Smedley  (1975) 
profile  unavailable,  R.  >,20 
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HEIGHr 


TABLE  III  (continued) 


C2 


(profile) 


Period  6 

4.19 

7,64 

13.89 

Period  7 

4.19 
6.6 
13.89 

Period  8 

4.19 
6.6 
13.89 

Period  9 

4.19 
6,6 
13.89 


28  March  74 


.060 
.030 
.024 


,014 

.0046 

* 


18  September  74 


,00030 
,00016 
.00014 


,0016 
,0012 
* 


19  September  74 


.00050 

.0025 

,00069 


.0063 
,0031 
.0015 


19  September  74 


.00062 

.0019 

.0011 


,0081 
.0036 
,0016 


2059-2117L 

5.9 
3.0 
2.3 

1643-1701L 

.029 
.015 
.013 

1735-1753L 

.048 

.24 

.07 

1900-1918L 

,059 

.18 

.11 


.09 
.19 
.52 


.12 
.27 
,68 


.11 
.23 

,64 


.04 
.08 
,09 
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Figure  17  includes  composites  of  height  versus  CJ;  for 
all  periods  for  both  the  stable  and  unstable  cases.   For 
these  composites,  obviously  unrepresentative  data  were  not 
plotted.   In  the  stable  case,  agreement  with  an  expected 
slope  greater  than  -3/2  is  possible.   The  lowest  level  shows 
the  greatest  range  in  relation  to  the  slope  line.   This  is 
attributed  to  a  comparatively  weak  amplitude  spectra  in 
many  cases.   Actually,  this  level  should  have  the  largest 
amplitude  spectra  in  comparison  to  the  other  three  levels. 
There  is  good  agreement  with  the  expected  -3/4  slope  for  the 
unstable  case. 
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VI ,   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

In  general,  the  height  variations  of  e ,  u+  and  C* 
observed  in  this  study  agree  with  the  overland  expressions 
developed  by  Wyngaard,   The  agreement  with  the  constant  flux 
assumption  is  an  important  result  of  this  study.   Most  devia- 
tions from  expected  results  were  able  to  be  attributed  to 
apparent  anomalies  in  spectral  results.   Correlation  of  C^, 
results  with  Richardson  number  derived  C™  values  was  not  good 
However,  since  spectral  results  were  generally  consistent 
between  levels,  the  spectral  results  are  considered  to  be 
more  reliable. 

It  has  been  shown  that  the  analog  spectral  analysis 
performed  in  this  study  produced  information  which  is  as 
reliable  as  that  from  digital  methods.   Application  of  the 
analog  method  on  the  large  amount  of  available  data  should 
furnish  information  as  to  the  reasons  for  identified 
departures  from  expected  results, 

B .  RECOMMEND AT I ONS 

Future  experiments  would  be  improved  by  relocating  level 
3  further  up  the  mast,  say  to  10  meters,  to  obtain  a  greater 
separation  between  levels  2  and  3. 

Improved  strip  chart  monitoring  of  u'  and  T1  signals 
during  experiments  is  needed.   Improved  monitoring  of  these 
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records  during  the  experiment  could  reveal  anomalous  signal 
levels  and  corrections  during  the  experiment  could  be 
accomplished.   In  general,  a  decrease  with  height   of  both 
u*  and  T'  signals  is  expected. 

More  frequent  in-situ  calibrations  of  the  hot  wire 
anemometers  should  be  made  to  increase  the  accuracy  of  results, 
especially  for  experiments  of  long  duration.   Frequent  resis- 
tance measurements,  on  the  platinum  wires  during  the  experi- 
ments is  also  necessary  for  this  purpose. 

Finally,  it  is  recommended  that  the  spectrum  analyzer's 
"vernier"  input  amplitude  control  be  tested  on  a  few  poor 
periods  to  determine  if  the  analyzer  sensitivity  could  be 
part  of  the  problem  in  the  lack  of  consistency  between  some 
levels.   Readjustment  for  each  level  and  also  obtaining  a 
spectral  amplitude  calibration  plot,  would  result  in  maximum 
sensitivity  and  minimum  distortion  of  the  spectral  output, 
A.C.  harmonics  from  the  ships  power  supply  should  also  be 
filtered  out  to  prevent . spectral  distortion. 
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APPENDIX  A 

A.   SPECTRUM  ANALYZERS-AVERAGER 

1 .  General 

The  following  discussion  is  for  further  amplification 
of  the  spectral  analysis  procedures  used  in  reducing  all  data 
for  this  study  and  will  serve  to  aid  those  who  will  continue 
to  study  turbulence  data  in  the  same  manner. 

2.  Equipment 

All  spectrum  analysis  was  performed  by  using  a 
Federal  Scientific  Model  UA-500-1  Spectrum  Analyzer-Averager, 
which  is  a  time-compression  spectrum  analyzer,  using  a  digital 
speed  up  loop . 

In  time-compression  spectrum  analyzers,  the  input 
data  is  digitized  and  stored  in  a  memory.   While  the  memory 
is  being  updated  with  new  data,  its  contents  are  read  out 
repeatedly  at  a  very  high  rate,   The  speeded-up  high  frequency 
replica  of  the  analog  input  signal  is  spectrum-analyzed,  using 
a  wide  filter  which  is  being  stepped  rapidly  in  real  time. 

A  recorder  analog  input  signal  is  converted  to  digital 
after  passing  through  an  anti-aliasing  low-pass  filter,   The 
sampling  frequency  is  controlled  by  the  analyzer  frequency 
operating  range.   The  digital  signal  is  stored  in  a  1500  word 
digital  memory,  circulates  in  the  memory,  and  can  be  non- 
destructively  read  out  at  an  accelerated  rate.   This  memory 
output  is  then  put  through  a  digital  to  analog  converter. 
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The  converter  output  is  a  time^compressed,  accelerated  form 
of  the  original  analog  input  signal,   It  passes  through  a 
heterodyner  and  is  applied  to  a  stable  crystal  filter,  causing 
the  filter  to  observe  successive  frequency-resolution  cells, 
A  weighting  function  modifies  the  analyzer's  frequency 
selectivity.   Settling  time  of  the  crystal  filter  is  short 
enough  to  obtain  one  spectral  point  each  time  the  memory  does 
one  recirculation.  Five  hundred  spectral  points  can  be  obtained 
in  50  milliseconds.   Averaging  of  the  spectral  wave  form  is 
then  accomplished  with  a  digital  ensemble  averager  to  obtain 
statistical  stability  in  mean  spectral  characteristics.   The 
final  result  in  storage  is  a  real-time  power-spectral- 
density  (PSD)  representation  available  at  output  terminals. 
The  averager  has  two  500  word  memories. 

A  Hewlett-Packard  X^Y  Plotter  was  used  to  obtain  a  hard 
copy  of  results  for  each  sampling  period. 

3,   Operating  Procedures 
a.   Analysis  Range 

The  range  considered  best  for  this  study  was  200 
Hz  for  all  data  analyzed.   This  resulted  in  a  sample  interval 
of  2.5  seconds  (memory  time).   Since  the  highest  frequency 
of  interest  was  about  100  Hz,  the  Nyquist  sampling  theorem 
was  satisfied.   The  sampling  rate  for  the  above  analysis 
range  is  600  Hz,  with  a  filter  bandwidth  of  0.4  Hz.   100 
percent  of  the  input  data  is  processed  for  the  200  Hz 
analysis  range. 
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b.  Input  Mode 

Input  signal  coupling  was  used  in  the  AC  mode  to 
analyze  only  the  AC  component  of  the  input  signal.   Hence, 
the  analysis  was  not  subject  to  DC  offsets.   This  mode  has 
a  low-frequency  3-dB  cutoff  at  less  than  3  Hz,   A  presampling 
filter  provides  at  least  50  dB  of  attenuation  per  octave 
fall  off  (i.e.  50  dB  at  400  Hz).   This  filter   and  a  sampling 
rate  of  three  times  the  analysis  range   eliminate  aliasing, 

c.  Input  Amplitude 

The  input  amplitude  was  adjusted  for  varying 
values  of  the  RMS  input  voltage  to  attenuate  the  level  of 
the  input  signal.   The  Vernier  Cal .  control  was  left  in  the 
detent  position  so  that  the  input  amplitude  switch  was 
functional  from  0-dB  through  -40-dB,  in  10-dB  steps  only. 
An  RMS  voltmeter  was  used  to  monitor  sample  portions  of  the 
input  analog  signal  to  determine  a  volts  RMS  full  scale. 
The  input  amplitude  was  then  adjusted  for  a  volt  setting  so 
that  full  scale  indications  didn't  exceed  this  volt  setting. 
Periodic  noise  spikes  were  not  considered.   Overflow  indica- 
tions were  not  normally  excessive  if  this  procedure  was 
followed.   This  indicator  may  flash  a  few  times  per  second 
without  danger  of  significant  signal  distortion. 

d.  Transient  Capture  Controls 

Depression  of  the  manual  button,  with  the  hold 
light  off,  enabled  the  memory  to  accept  new  data  continuously 
The  two  other  available  modes  were  not  used  in  this  analysis. 
Then  the  trigger  polarity  and  sensitivity  controls  and  the 
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transient  auto  delay  sv/itches  are  not  functional,   The 
external  sampling  switch  must  also  be  off, 

e.  Output  Spectrum  Gain 

Pushbutton  operation  introduces  a  gain  (0  dB, 
+10  dB,  +  20  dB)  to  the  analyzer  output  prior  to  averaging  of 
spectra.   Operation,  for  this  study,  was  in  the  +10  dB  (gain 
x  3.2)  mode,  due  to  the  characteristics  of  the  data  signals. 
The  output  gain  overload  light  indicates  when  spectral  peaks 
are  above  full  scale  and  are  clipped,   This  occurred 
infrequently . 

f.  Amplitude  and  Frequency  LIN/LOG  Switches 
Logarithmic  scale  for  both  amplitude  and 

frequency  were  used  for  all  spectra. 

g.  Display  Sweep  Controls 

The  FAST  OSCILL  mode  must  be  used  for  the  averaging 
of  a  signal  and  the  SLOW  PLOT  mode  should  be  selected  for 
plotting  the  memory  contents  of  the  X-Y  plotter.   The  START, 
STOP,  RESET  pushbuttons  .control  the  averaging  or  plot  mode 
elected.   In  slow  plot  an  automatic  rate  control  is  used  for 
large  signal  excursions  to  avoid  plotter  slew-rate  problems, 
h.   Averager 

(1)  Active  Memory  Control,   This  switch  selects 
that  memory  where  averaging  takes  place  and  whose  contents 
are  displayed  or  plotted.   The  ERASE  pushbutton  will  erase 
the  selected  memory  if  FAST  OSCILL  is  selected. 

(2)  Averaging  Method  Controls.   In  the  SUM  mode, 
each  non-redundant  analyzer  spectrum  is  added  to  previous 
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spectra  stored  in  the  active  memory,   The  INSTANT  SPECT 
control  can't  be  depressed  when  an  averaging  process  is  to 
be  started. 

(3)   Average  Stop  Mode  Controls,   The  manual 
mode  was  used  in  this  study  to  stop  each  analysis  at  18 
minutes  with  the  STOP  control.   The  Number  of  Spectra-N 
switch  then  functioned  as  a  gain  control  if  set  at  N  less 
than  256.   This  proved  to  be  a  problem,  in  that  the  analyzer 
looks  at  this  switch  to  indicate  how  many  samples  are  coming, 
If  the  number  N  selected  and  the  number  N  of  samples  do  not 
coincide  when  analysis  stops,  then  the  spectra  is  not 
complete.   There  is  no  problem  in  the  Stop-At-N  mode. 
i.   Calibration 

All  other  operating  procedures  and  the  calibration 
procedures  of  Analyzer-Averager  and  displays  can  be  obtained 
from  the  manufacturer's  instruction  manual, 
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APPENDIX  B 

A.   SPECTRAL  X-Y  PLOT  SCALING  PROCEDURES 

To  obtain  a  power  spectral  density  level,  corresponding 
to  an  RMS  voltage  input,  a  calibrated  scale  overlay  had  to 
be  constructed  for  use  in  spectral  interpretation. 

Amplitude  scaling  was  accomplished  by  an  internally 
generated  signal  available  in  the  AMPL,  CAL  input  mode, 
When  O-dB  (,lv)  input  gain  was  selected,  a  cal  signal  with 
amplitude  approximately  equal  to  0,1  volts  RMS  appears  on 
the  display  as  a  spike  at  100  Hz  indicating  full  scale 
amplitude.   This  spike  was  then  plotted  in  the  X-Y  plotter, 
along  with  bench  marks,  using  a  spectral  gain  of  OdB  (xl) 
setting.   Successively,  the  input  gain  was  stepped  to 
attenuate  the  amplitude  of  the  spike  by  10  dB,  and  a  plot  of 
the  height  of  each  spike  was  obtained  on  the  X-Y  plot.   This 
plot  now  represented  a  graduation  of  RMS  input  from  a  minimum 
of  0,001  volts  to  0.1  volts  maximum.   These  voltages  were 
then  converted  to  LOG  10  and  a  graduated  scale  was  constructed 
to  that  the  LOG  10  of  volts  RMS  could  be  interpolated  from  a 
spectral  plot.   The  amplitude  scale  was  adjusted  for  each 
spectrum  as  a  function  of  input  attenuation  and  spectrum  gain. 
These  LOG  10  values  were  then  converted  to  PSD  levels  for  use 
in  calculating  epsilon  and  structure  parameter  values.   The 
Federal  Scientific  Manual's  equation  to  convert  to  PSD  is 
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PSD  level  (V2/Hz)  =  °-79  <ca1'  level  VRMS)2 

Filter  Bandwidth 

where  the  filter  bandwidth  ($  =  0.4)  is  a  function  of  the 
analysis  range  selected, 
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